CERN-PH-TH/2007-128 



Collinear improvement of the BFKL kernel in 

the electroproduction of two light vector 

mesons 



F. Caporale 1 , A. Papa 1 and A. Sabio Vera 2 

1 Dipartimento di Fisica, Universitd della Calabria 
and Istituto Nazionale di Fisica Nucleare, Gruppo collegato di Cosenza 
1-87036 Arcavacata di Rende, Cosenza, Italy 

2 Physics Department, Theory Division, CERN, 
CH-1211, Geneva 23, Switzerland 



February 5, 2008 



Abstract 

The use of the BFKL kernel improved by the inclusion of subleading terms generated by 
renormalization group (RG) analysis has been suggested to cure the instabilities in the 
behavior of the BFKL Green's function in the next-to-leading approximation (NLA). We 
test the performance of a RG-improved kernel in the determination of the amplitude of a 
physical process, the electroproduction of two light vector mesons, in the BFKL approach in 
the NLA. We find that a smooth behavior of the amplitude with the center-of-mass energy 
can be achieved, setting the renormalization and energy scales appearing in the subleading 
terms to values much closer to the kinematical scales of the process than in the approaches 
based on the unimproved kernel. 



1 Introduction 



It is known that hard processes in which the center-of-mass energy is much larger than 
all the other scales are the natural ground for the application of the BFKL approach [1]. 
This approach was originally developed in the leading logarithmic approximation (LLA), 
which means resummation of all terms of the form (a s ln(s)) n . In such an approximation the 
argument hr of the running coupling and the energy scale are not fixed. This motivated the 
extension of the approach to the next-to-leading logarithmic approximation (NLLA), which 
means resummation of all terms proportional to a s (a s m(s)) n . In both approximations 
the BFKL amplitude appears as a convolution of the Green's function of two interacting 
Reggeized gluons with the impact factors of the colliding particles (see, for example, Fig. 1). 
The Green's function, which carries the dependence on the center-of-mass energy, can be 
determined through the BFKL equation. The impact factors are process-dependent and 
describe the interaction between Reggeized gluons and scattering particles. 

The singlet kernel of the BFKL equation in the next-to-leading approximation (NLA) was 
obtained for the forward case in Ref. [2], completing the long program of calculation of the 
NLA corrections [3] (for a review, see Ref. [4]). In the non-forward case the ingredients for the 
NLA BFKL kernel have been known for a few years in the case of the color octet representa- 
tion in the t-channel [5]. This color representation is very important to check the consistency 
of the s-channel unitarity with the gluon Reggeization, i.e. for the "bootstrap" [6]. More 
recently, the last missing piece for the determination of the non-forward NLA BFKL kernel 
has been calculated in the singlet color representation, i.e. in the Pomeron channel, relevant 
for physical applications [7] . The singlet NLA BFKL kernel in the so-called "dipole form" is 
available now also in the coordinate representation [8], which allows the study of its confor- 
mal properties and the comparison with the kernel of the Balitsky-Kovchegov [9] equation 
in the linear regime. So far, the color dipole kernel has been calculated in the NLA only for 
the quark part [10] and agrees with the dipole form of the quark part of the NLA BFKL 
kernel. 

In this paper we will focus on the BFKL approach in the NLA and in the case of forward 
scattering. It is well known that the NLA corrections to the Green's function turn out 
to be large, this being a signal of the poor convergence of the BFKL series. In order to 
"cure" the resulting instability, more convergent kernels have been introduced, including 
terms generated by renormalization group (RG), or collinear, analysis [11]. They are based 
on the upshift method [11], with uj being the variable Mellin-conjugated to the squared 
center-of-mass energy s. The main effect of this method is that the scale-invariant part of 
the kernel eigenvalues carries a dependence on the Mellin variable u, in such a way that the 
position of the singularities of the Green's function in the cu-plane becomes the solution of an 
implicit equation in uj. Many other studies have been performed, either based on this kind 
of improved kernels [12] or analyzing different aspects of the kernel NLA and alternative 
approaches [13]. The effects of these collinear corrections in exclusive observables have been 
investigated in Ref. [14], with a posteriori confirmation in Ref. [15]. 

In Ref. [16] the original approach of Ref. [11] was revisited and an approximation to 
the original upshift was performed, leading to an explicit expression for the RG-improved 
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NLA kernel. It was shown that this improved kernel leads to a NLA BFKL Green's function 
exempt of instabilities. Since the effect of the RG-improvement is to modify the BFKL kernel 
by the inclusion of terms beyond the NLA, one is led to conclude that RG-generated terms, 
although formally subleading, play an important numerical role in practical applications. 

It is very interesting to test the RG-improvement of the kernel in the calculation of a full 
physical amplitude, rather than just considering its effect on the BFKL Green's function, 
and to compare it with other approaches. A test-field for this comparison can be provided by 
the physical process 7*7* — > VV, where 7* represents a virtual photon and V a light neutral 
vector meson (p°,u,(/>). The amplitude of this reaction has been calculated in Ref. [17] 
through the convolution of the (unimproved) BFKL Green's function with the 7* — > V 
impact factors, calculated in Ref. [18] 1 . In the case of equal photon virtualities, the so- 
called "pure" BFKL regime, a numerical calculation has shown that NLA corrections are 
large and of opposite sign with respect to the leading order and are dominated, at the lower 
energies, by the NLA corrections from the impact factors. Nonetheless, an amplitude for 
this process with a smooth behavior in s could be achieved by "optimizing" the choice of 
the energy scale s and of the renormalization scale /i R , which appear in the subleading 
terms. Later on it has been found that the result is rather stable under change of the 
method of optimization of the perturbative series and of the representation adopted for the 
amplitude [19]. 

The striking feature of these investigations was that in all cases the optimal values of the 
two energy parameters turned out to be quite far from the kinematical scales of the reaction. 
For example, the optimal value of the renormalization scale hr turned out to be typically as 
large as ~ 10Q, Q 2 being the virtuality of the colliding photons. The proposed explanation 
for these "unnatural" values was that they mimic the unknown next-to-NLA corrections, 
which should be large and of opposite sign respect to the NLA in order to preserve the 
renorm- and energy scale invariance of the exact amplitude. If this explanation is correct and 
if the RG-improvement of the kernel catches the essential dynamics from subleading orders, 
then, by repeating the numerical determination of the 7*7* — > VV amplitude with the use of 
an RG-improved kernel, one should get more "natural" values for the optimal choices of the 
energy scales and, of course, results consistent with the previous determinations. In this work 
we address this question by calculating the NLA amplitude of the 7*7* — > VV process in the 
BFKL approach with the RG-improved kernel of Ref. [16], which can be straightforwardly 
implemented in the numerical set up of Refs. [17, 19]. 

The paper is organized as follows: in the next Section we repeat the steps of Refs. [17, 
19] to build up the NLA amplitude in two representations, series and "exponentiated", 
which implement the RG-improved kernel of Ref. [16]; in Section 3 we numerically evaluate 
the amplitude, considering both the cases of colliding photons with the same virtualities 
and with strongly ordered virtualities. We stress that in Refs. [17, 19] only the case of 
equal photons' virtualities was considered; attempts to determine the amplitude for strongly 
ordered virtualities were unsuccessful, due to the large instabilities met in the numerical 
analysis [20]. We expect that the RG-improvement should be even more effective in the 
latter case, since it was conceived to work in a kinematics with strong asymmetry in the 

1 This amplitude has been considered also in [22, 23, 24]. 
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transverse momentum plane [11]. 



2 The NLA amplitude with the RG-improved Green's 
function 2 

We consider the production of two light vector mesons (V = p°, u>, 0) in the collision of two 
virtual photons, 

7*(p) 7*(p') - V( Pl ) V(p 2 ) . (1) 

Here, p\ and p 2 are taken as Sudakov vectors satisfying p\ = p\ = and 2(p 1 p 2 ) = s; the 
virtual photon momenta are instead 

p = api p 2 , p =ap 2 r pi , (2) 

as a's 

so that the photon virtualities turn to be p 2 = —Q\ and {p') 2 = —Q\. We consider the 
kinematics when 

s > Q\ 2 > A 2 QCD , (3) 

and 

a = 1 + 91 + 0( s ~ 2 ) , = 1 + 5! + 0(s- 2 ) . (4) 

In this case the vector mesons are produced by longitudinally polarized photons in the 
longitudinally polarized state [18]. Other helicity amplitudes are power suppressed, with a 
suppression factor ~ mv/Qi, 2 - We will discuss here the amplitude of the forward scattering, 
i.e. when the transverse momenta of the produced V mesons are zero or when the variable 
t = (pi — p) 2 takes its maximal value t = —Q\Q\js + 0(s~ 2 ). 

The forward amplitude in the BFKL approach may be presented as follows 

d 2 ~* d 2 ~* <5+*co ^ w 

^^^/f^'^/f^^ 881 / 2^ (£) (5) 

<5— ioo 

This representation for the amplitude is valid with NLA accuracy. 

In Eq. (5), $i(gi,so) and § 2 {— q 2 ,so) are the impact factors describing the transitions 
7*(p) — > V(pi) and 7*(p') — > V(p2), respectively. The Green's function in (5) is determined 
by the BFKL equation 

<5 2 (<fi - 92) =wG w (gi,g 2 ) - / d 2 qK(q 1} q)G u (q,q 2 ) , (6) 



2 This Section follows closely Section 2 of the Refs. [17, 19], the only difference being the use of a modified 
BFKL kernel. The reader already familiar with the notation and the previous papers may prefer to go 
straight to the main formulas: Eq. (36) for the "exponentiated" representation, Eq. (37) for the "series" 
representation of the amplitude and Eq. (22) for the extra-term in the BFKL kernel eigenvalue. 
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Figure 1: Schematic representation of the amplitude for the 7*(p)7*(p') — > V(pi) V(p 2 ) 
forward scattering. 

where K(qi,q 2 ) is the BFKL kernel. It is convenient to work in the transverse momen- 
tum representation, where "transverse" refers to the plane orthogonal to the vector mesons 
momenta. In this representation, defined by 

= ( 7 ) 
(?i|g 2 > = 5 (2) (<?i " <f 2 ) , = (A\k)(k\B) = J d 2 kA(k)B(k) , (8) 

the kernel of the operator K is 

K{q 2 ,q l ) = (q 2 \K\q l ) (9) 
and the equation for the Green's function reads 

i = (u - k)G u , (io) 

its solution being 

G u = (u-K)- 1 . (11) 
To clearly indicate the RG-improved pieces of the kernel, we decompose K as 

K = a s K° + a 2 s K l + K RG , (12) 

where 

a s = (13) 

7T 

and N c is the number of colors. In Eq. (12) K° is the BFKL kernel in the LLA, K 1 is the 
NLA correction and Krq includes the RG-generated terms, which are 0{a\). The impact 
factors are also presented as an expansion in a s 



^iM = oc s D h2 \Ci%q 2 ) + a s Ci%q 2 )\ , D lfl = -^%±y/N* -1 , (14) 



Aixe q f v 
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where fy is the meson dimensional coupling constant (f p ~ 200 MeV) and e q should be 
replaced by e/y/2, e/(3y/2) and — e/3 for the case of p°, to and meson production, respec- 
tively. 

In the collinear factorization approach the meson transition impact factor is given as 
a convolution of the hard scattering amplitude for the production of a collinear quark- 
antiquark pair with the meson distribution amplitude (DA). The integration variable in this 
convolution is the fraction z of the meson momentum carried by the quark (z = 1 — z is the 
momentum fraction carried by the antiquark): 

^'hwrkr^- (15) 

o ' 

The NLA correction to the hard scattering amplitude, for a photon with virtuality equal to 
Q 2 , is defined as follows 

c,1,(?2 > = wj a dz i4te? 1t(:) + T<1 - z)]Uz) • (16) 

with t(z) given in the Eq. (75) of Ref. [18]. C^itf 2 ) are given by the previous expression 
with Q 2 replaced everywhere in the integrand by Q\ and Q 2 ,, respectively. We will use the 
DA in the asymptotic form 4>^ s (z) = 6z(l — z). 

To determine the amplitude with NLA accuracy we need an approximate solution of 
Eq. (11). With the required accuracy this solution is 



G u = (u - asK )- 1 + (u- ask*)- 1 (a^K 1 + K RG ) {uo - a s k Q y l + O U^K 1 ) 



• (17) 



Differently from Refs. [17, 19], where K R g wa s absent, this Green's function includes effects 
which are beyond the NLA. The basis of eigenfunctions of the LLA kernel, 

K°\v) = X {u)\u) , x iy) = 2^(1) - V> (I + w) - V> (I - , (is) 

is given by the following set of functions: 

(9>> = ^ (« 2 ) i "" i ' (19 » 

for which the orthonormality condition takes the form 

<-»=/0(5 2 )"'-'"'- 1 = ^-.'). (20) 



2tt 2 

The action of the modified BFKL kernel on these functions may be expressed as follows: 



+ ^r)^-xW) \y) + Xrg{v)\v) , (21) 



where the first term represents the action of LLA kernel, the second and the third ones stand 
for the diagonal and the non-diagonal parts of the NLA BFKL kernel [17] and 



XRG {u) = 2Ke £ 

lm=0 



E 



-l) n (2n)! 



/ — . -2\ n + 1 



a. 



k^o 2 n n\(n + 1)! (1/2 + iv + m - b a 8 ) 
_ ( a b 



2n+l 



(22) 



+ 



1/2 + iv + m s \l/2 + iv + m (1/2 + iv + m) 2 2(1/2 + iv + m) 3 1 

is the solution of the cu-shift equation obtained in [16], with 

5 fio 13 71/ 55 1 /?o ?V 11 

a = 12 iVc ~ 36 iVf ~ 36' = ~ 8 W c ~ 6iVf ~ 12' 

The function * s conveniently represented in the form 



(23) 



= -g|r (x» - yxW " *xV)) + XW , 



(24) 



where 



^- i x(z,)-6C(3)-x»- 7 



+ 



1 

<f>(v) = 2 J dx 



n 2 smh(nv) 
2 z/ cosh 2 (7rz/) 

cos(i/ln(x)) 



3+1 + -^ 



C0Sh(7Tf 

71/ \ ll + 12z/ 2 ^ 



a" c 3 / 16(1 + z/ 2 ) 



+ 4 0(z/) 



(l + a;)y / a : 6 



7T 



— - U 2 (x) 



Li 2 (x) = - 1 



ln(l - t) 



(25) 
(26) 



Here and below x'M = d(x{v))/dv and x"^) = d 2 {x{v)) / 'd 2 v. 

The 1 1/) representations for the impact factors are given by the following expressions: 
C[°\q 2 ) + ?° 



dv' Cl {v'){v'\q) , 



Cf(9 2 ) +? ° 



= | rfz/c 2 (z/)(g» , (27) 



(? 2 ) 



(? 2 ) 



2\-«f- 



(28) 



and by similar equations for \v) and c^ 1 ^) from the NLA corrections to the impact 
factors, C{ 1] (q 2 ) and C?\q 2 ). 

Following Ref. [17], we obtain the amplitude as a spectral decomposition on the basis of 
eigenfunctions of the LLA BFKL kernel: 

+oo 



D X D 2 (2nf 



s 



1 dv[ — \ a 2 s {^ R ) Cl {v)c 2 {v) ll + a s {^i R ) ' 



+ « 2 (^i?) ln ( — 



ci(z/) c 2 (z/) 

10 rfln(^) 

lc2Mj +21n(^) 



+ In ( - ) X rg(v) } ■ 



(29) 



We find that 



ci(v)c 2 {v) = 



(Ql 



±iu- 



r 2 [|±H 6 



7T 



y/2 r[3±2iz/]cosh(7rz/) ' 
1 (QlY 97r 3 (l + 4i/ 2 )sinh(7Ti/) 



i ^ n UoV = 2 
dv 



Q1Q2 \Ql) 32 v (1 + J/ 2 ) cosh 3 (7rz/) ' 
V>(3 + 2ii/) + ^(3 - 2ii/) - V ( I + iv) ~ i> ( I ~ iv) ~ In (Q1Q2) 



(30) 
(31) 
(32) 



It can be useful to separate from the NLA correction to the impact factor the terms 
containing the dependence on s and fa, 
1 



dz 



q 



q 2 + zzQ 2 



0|| U 



(33) 



x 



In 



Sq_ 

,Q 2 



In 



(a + zz) 4 " 
a 2 z 2 z 2 



+ 



_Po_ 

4:N r 



In 



Q 2 ) ' 3 



+ - - ln(a) + 



Accordingly, one can write 

c§(v) = + (34) 

where c^(u) are the contributions from the terms isolated in the previous equation and 
Ci^iv) represent the rest. In Ref. [17] it was found that 



ci>) + 4>) 
d{v) c 2 {y) 



In 



so 



Q1Q2 



2N r 



In 



/4 
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Q1Q2J + 3 



+ ^(3 + 2ii/) + ^(3 - 2iv) -ipl-z + iv) -^(--iv 



(35) 



One can construct infinitely many representations of the amplitude, all of them equivalent 
within NLA accuracy. A particular one, motivated in Ref. [19], is to exponentiate all the 
scale-invariant part of the NLA kernel, obtaining 



Xm s {A) s 
D X D 2 = (2^)2 

X 1 + OC a (fJ, R ) 



+oc 



dv[ — 

.so 



a s 2 (/ifi)ci(^)c 2 (^) 



ci(z/) c 2 (z/) 



+ 2 ln(/i 



(36) 

Another possible representation of the amplitude, in some sense closer to the original idea 
of the BFKL approach, is the "series" representation, which reads 

Q1Q2 2m s A 1 



D X D 2 s 



x 



f / s \ 00 / s \ n+l 

\b + a In — + W a n In — (37) 

I Vsq/ n= i Vso/ 



Vs 



+ 6 n In — + d n (s ,ijL R ) In [ — 



\s 



n-l 
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where the coefficients 



+00 

J dvci(v)c 2 (v) 



Q1Q2 

—00 

are determined by the kernel and the impact factors in LLA and 



'1Q2 J n\ 



Q 



arise from the collinear improvement. The coefficients 



d n — n In 



so 



Q1Q2 



+ 



f n+ l)^ziin 



fj? R \ n(ra-l) 



Q1Q2 



(38) 



(39) 



Q1Q2 



dv (n + I) f {y)ci{y)c 2 (v 



Q1Q2 



+00 



dv ci(v)c 2 {v) 



(n-l)\ 



n- 1)! 

c«(^ , %\v) 



+ 



ci(u) c 2 (v) 



+ (n- 1) 



xiy) 
xiy) 



(40) 



are determined by the NLA corrections to the kernel and to the impact factors. Here, c^ 2 (v) 
represent the contribution without the terms depending on s and (3 , and 



f(v) = - + ^(3 + 2iv) + ^(3 - 2ii/) -^(2+^-^(2- • 



(41) 



We stress that the terms in the series representation (37) with the a n coefficients are beyond 
the NLA, since, as one can easily see from Eq. (22), xrg is 0(a^). 



3 Numerical results 



In this section we present some numerical results for the dependence in s of the BFKL 
amplitude calculated for the process under study, using both the "exponentiated" and the 
"series" representations derived in the previous Section. Following Ref. [17], we will adopt 
the principle of minimal sensitivity (PMS) [21] requiring, for each value of s, the minimal 
sensitivity of the predictions to the change of both the renormalization and the energy scale, 
Hr and So- In previous studies, where the unimproved kernel was used, the optimal choices 
for jiR and So turned out to be very far from the kinematical scales of the process. Our 
aim is to see if and to what extent the inclusion of a collinear improvement leads to more 
"natural" values for the optimal scales. This would demonstrate that the RG-generated 
terms reproduce the essential subleading dynamics, thus stabilizing the perturbative series. 
In the following analysis we use the two-loop running coupling corresponding to the value 
a s (M z ) = 0.12. 
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Figure 2: Tm s (A)Q 2 / (s DiD 2 ) as a function of Y at Q 2 =24 GeV 2 and rif = 5 in the 
"exponentiated" representation with and without collinear improvement of the kernel; in 
both cases the PMS optimization method has been used. 

3.1 Symmetric kinematics 

We consider here the Qi = Q 2 = Q kinematics, i.e. the "pure" BFKL regime, with 
Q 2 =24 GeV 2 and rif = 5. We start with the "exponentiated" representation, given in 
Eq. (36) and set ln(s/s ) = Y - Y , where Y = ln(s/Q 2 ) and Y = \n(s /Q 2 ). We have 
looked for the optimal value for the scales /ir and Y . In practice, for each fixed value of Y 
we have determined the optimal choice of these parameters for which the amplitude is the 
least sensitive to their variation. We have found that the amplitude is always quite stable 
under variation of both scales and exhibits generally only one stationary point (local maxi- 
mum). We choose as optimal values of the parameters those corresponding to this stationary 
point. 

The optimal values turned out to be typically hr ~ 3Q and Y ~ 2. In comparison 
with Ref. [17], where the optimal choice was typically /j,r ~ 10Q, we can see that there is a 
remarkable move towards "naturalness" . The fact that the inclusion of the RG-terms affects 
the optimal choice of hr more strongly than of Y is not surprising, since the added terms 
depend on /ir and not on Y . In Fig. 2 we show the result for the (imaginary part of the) 
"improved" amplitude compared with the result obtained in Ref. [19]. The curves are in 
good agreement at the lower energies, the deviation increasing for large values of Y. This 
is consistent with having a larger asymptotic intercept when the collinear improvements are 
taken into account. We have to remember, however, that the applicability domain of the 
BFKL approach is determined by the condition 6i s (hr)Y ~ 1, that for our typical optimal 
value of /iR and for Q 2 =2A GeV 2 means Y ~ 6. Around this value the discrepancy is not so 
pronounced. 
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Figure 3: lm s (A)Q 2 / (s DiD 2 ) as a function of Y at Q 2 =2A GeV 2 and rif = 5 in the "series" 
representation with and without collinear improvement of the kernel; in both cases the PMS 
optimization method has been used. 

The next analysis has been done using the "series" representation of the amplitude, given 
in Eq. (37). In this case we have also observed a smooth dependence of the amplitude on 
the two energy parameters. The optimal values for Y and turned out to be quite similar 
to those obtained for the "exponentiated" representation, ~ 3Q and Y ~ 3. In Fig. 3 
we show the behavior in Y of the "series" amplitude, compared with the determination of 
Ref. [17]. The situation is similar to Fig. 2, but the deviation between the curves appears 
to be more marked here. It is important to observe that the curves for the "exponentiated" 
and "series" representations of the amplitude as functions of Y with collinear improvement 
(see Figs. 2 and 3) fall almost on top of each other, while in the determination without the 
collinear improvement there was a discrepancy, more pronounced at higher energies [19]. 
This is a further indication of a better stability, induced by the collinear improvement. 

In order to make visible the effect of the collinear improvement in the "series" represen- 
tation we list the first few coefficients (see Eq. (37)) b n , d n , coming from the unimproved 
BFKL kernel and impact factors (in LLA e NLA respectively), and a n , coming from the 
RG-resummed terms. Using the optimal scales chosen with the PMS method we obtain 
(Q 2 =24 GeV 2 , n f = 5, Y = 3, n R = 3Q) 

b = 17.0664 6i = 34.5920 b 2 = 40.7609 b 3 = 33.0618 b A = 20.7467 

di = 0.674275 d 2 = -1.73171 d 3 = -7.46518 d 4 = -15.927 (42) 
ai = 5.52728 a 2 = 7.30295 a 3 = 6.42149 a 4 = 4.24011 . 

We can see that the a n coefficients are of the opposite sigh respect to the d n) so "curing" 
the bad behavior of the BFKL series. Even if the values of the a n coefficients go down 
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Figure 4: Xm s (A)QiQ2/{s D1D2) as a function of Y for photons with strongly ordered 
virtualities (Q2/Q1 = 6 and Q2/Q1 = 96, with QiQ2=24 GeV 2 ), in comparison with the 
case of photons with equal virtualities {Q\ = Q\=2A GeV 2 ). All curves have been obtained 
using the "exponentiated" representation with the collinearly improved kernel. 

with n, they appear in Eq. (37) with two more powers of the energy logarithm than the d n 
coefficients, so that their effect is not limited to low energies. 

3.2 Asymmetric kinematics 

When the virtualities of the photons are strongly ordered, we enter the "DGLAP" regime, 
where collinear effects should come heavily into the game. In this regime, previous attempts 
to numerically determine the amplitude using unimproved kernels were unsuccessful due to 
severe instabilities [20] . We have found here that these instabilities disappear if, instead, the 
RG-improved kernel is used. 

In the numerical analysis to follow, we consider two choices for the virtualities of the pho- 
tons, Qi=2 GeV, Q 2 =12 GeV and Qi=0.5 GeV, Q 2 =48 GeV, so that Q X Q 2 = Q 2 =24 GeV 2 
in both cases, and used the "exponentiated" representation. We define Y = ln(s/QiQ 2 ) and 
Y = \n(s /Q 1 Q 2 ). 

For the first choice of virtualities, we find that for each Y value the amplitude is still quite 
stable under variation of the energy parameters and the optimal values are /ir ~ AVQ1Q2 
and Y ~ 2, almost independently of Y. The same holds for the second choice of virtualities, 
with the only difference that now the optimal values depend strongly on Y. As an example, 
for Y = 6, when a s ([iR)Y ~ 1, the optimal fiR is ~ 3^/QiQ2, but Y =7. This large value for 
Yq should not be surprising: if we use Q\ as normalization scale in Yq instead of Q1Q2, the 
optimal value lowers down ~2.5, which looks more "natural". 
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In Fig. 4 we plot the amplitude for the two choices of photons' virtualities we have 
considered, together with the amplitude for Qi = Q 2 = \p2A GeV. The amplitude becomes 
smaller and smaller when Q2/Q1 increases, as it must be expected due to the presence of the 
factor cos(z/log(<52/<3i)) in the integration over v. We stress again that, if the RG-generated 
terms are removed, it is impossible even to draw the curves in Fig. 4 with Q 2 ^ Qi- 



4 Conclusions 

We have applied a RG-improved kernel to determine the amplitude for the forward transition 
from two virtual photons to two light vector mesons in the Regge limit of QCD with next- 
to- leading order accuracy. The result obtained is independent on the energy scale s , and 
on the renormalization scale /ir within the next-to-leading approximation. 

Using two different representations of the amplitude, which include the dependence on 
the energy scale and on the renormalization scale at subleading level, we have performed a 
numerical analysis both in the kinematics of equal and strongly ordered photons' virtualities. 

An optimization procedure, based on the principle of minimal sensitivity, has led to 
results stable in the considered energy interval, which allow to predict the energy behavior 
of the forward amplitude. The important finding is that the optimal choices of sq and hr are 
much closer to the kinematical scales of the problem than in previous determinations based 
on unimproved kernels. This effect is very marked for fiR, as it must be expected, since the 
extra-terms depend on fi R and not on s . This leads us to conclude that the extra-terms 
in the BFKL kernel coming from collinear improvement, which are subleading to the NLA, 
catch an important fraction of the dynamics at higher orders. 

Moreover, the use of the improved kernel has allowed to obtain the energy behavior of 
the forward amplitude in the case of strongly ordered photons' virtualities, which turned out 
to be unaccessible to previous attempts using unimproved kernels. 
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